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Abstract 

Using molecular dynamics, we study the nucleation and stability of bulk nanobubble clusters. We study 
the formation, growth, and final size of bulk nanobubbles. We find that, as long as the bubble-bubble 
interspacing is small enough, bulk nanobubbles are stable against dissolution. Simple diffusion calculations 
provide an excellent match with the simulation results, giving insight into the reason for the stability: 
nanobubbles in a cluster of bulk nanobubbles "protect" each other from diffusion by a shielding effect. 



1 Introduction 

Gas bubbles arc ubiquitous in nature, industry and 
daily life. They are found in streams of water, manu- 
facturing processes of many types of materials, and, 
of course, when we enjoy a carbonated drink. Even 
though gas bubbles are commonly present in many 
of the liquids we deal with on a daily basis, bubbles 
are, in fact, usually unstable against dissolution in 
the medium that surrounds them [1]. The dissolution 
rate increases as the bubble becomes smaller because 
of the increased (Laplace) pressure Ap = 27/i? in- 
side the bubble, where 7 is the interfacial tension of 
the bubble wall and R the bubble radius. The conse- 
quence is that nanoscopic bubbles cannot survive for 
more than a few microseconds. 

In contrast to this expectation, surprisingly, exper- 
iments by Ohgaki et al. [2] have shown that stable 
bulk nanobubbles do exist. In these experiments the 
bubbles were observed to be packed closely together 
(the distance between neighbouring bubbles was mea- 
sured to be less than lOi?), suggesting that a shield- 
ing mechanism between bubbles may act to keep the 
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bubbles from dissolving. In addition to this direct 
observation of bulk nanobubbles, their presence has 
also been indirectly measured in experiments, using 
dynamic laser light scattering [3,4]. Although this 
technique cannot distinguish between nanobubbles 
and liquid density variations in the liquid caused by 
other sources (such as large organic molecules), the 
observed fluctuations disappear after degassing the 
liquid, indicating that the observed objects are in- 
deed bulk nanobubbles. 

In addition to these experiments, there are many 
publications where the presence of surface nanobub- 
bles are observed at liquid-solid interfaces. Generally, 
these surface nanobubbles are detected by Atomic 
Force Microscopy (AFM), and they can survive for 
days [5-7]. Similar to bulk nanobubbles, surface 
nanobubbles should dissolve within microseconds, in 
contrast to the AFM observations [5,6,8-17]. Vari- 
ous stabilization mechanisms have been proposed [5- 
7, 18, 19], and many of them invoke the direct bubble- 
wall interaction. This in particular holds for the 
dynamic equilibrium theory promoted by some of 
us [7, 19,20]. This stabilization mechanisms is there- 
fore not applicable to bulk nanobubbles: the symme- 
try breaking caused by the presence of the substrate 
in the case of surface nanobubbles does not exist for 
bulk nanobubbles. 
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On the other hand, different stabiUzation mecha- 
nisms may exist that could account for stable bulk 
nanobubblcs. Such a mechanism will be discussed in 
this paper: when a bulk nanobubble is surrounded by 
more nanobubbles, the diffusive outflux is 'shielded': 
a locally high concentration of dissolved gas in the 
water suppresses the diffusive outflux from the bub- 
ble. For this to happen a cluster of bubbles must exist 
where the spacing between bubbles is not too large. 
Indeed, the bulk nanobubbles reported by Ohgaki et 
al. [2] have a distance of lOR or less. 



In this paper, we will discuss molecular dynam- 
ics (MD) simulations of binary mixtures of simple 
(Lennard- Jones) fluids. One of the fluids is under 
the imposed conditions (T = 300K, p = 10^ Pa) 
in the liquid state, the other in the gaseous state. 
The simulations will be carried out in a simula- 
tion domain of which one dimension is very small 
{£ X i X d, d ^ £) , such that the simulations are quasi- 
2D, see also Fig. 1. For a full 3D case, the results 
will only differ quantitatively, but qualitatively they 
will be the same. Periodic boundary conditions are 
applied in all directions, such that we only have to 
simulate one single nanobubble which is then mir- 
rored. This infinite repetition of nanobubbles then 
represents an infinite (periodic) nanobubble cluster 
( see Fig. 1) in a closed system. The closed system 
means that the total amount of gas is conserved. In 
this work, we explore two box sizes: £ ~ 15 nm and 
£ = 30 nm, with d — 3.64 nm in both cases. 



The paper is organized as follows: In section 2 the 
numerical details of the simulations will be outlined, 
such as the parameters and algorithms used, as well 
as the initial conditions. Next, in section 3, the re- 
sults of the MD-simulations will be presented and 
discussed, and in section 4 we will compare these re- 
sults (particularly the equilibrium radius R^q of bulk 
nanobubbles) with continuum predictions. Finally, 
in section 5 we will discuss the stability of the entire 
cluster (as opposed to just single bubbles inside the 
cluster) . 




Figure 1: A cluster of nanobubbles on a rectangular 
grid. The dotted lines indicate the unit cell, which is 
rectangular and has sides with length I: the distance 
between two neighbouring bubbles. 



2 Numerical details 

2.1 Molecular Dynamics 

To simulate a cluster of bulk nanobubbles, we use 
Molecular Dynamics (MD) simulations of simple flu- 
ids. The atoms in the simple fluids interact with each 
other through the Lennard- Jones potential: 



ULjir) = 4e, 



r J \ r 



(1) 



Here, is the interaction strength between atom 
species i and j, and cTy the interaction radius 
between atoms species i and j. In our simu- 
lations we use two atom types: the first is in 
the liquid state under the conditions considered 
{p = 10^ Pa, T = 300K) and the second in 
the gas state. The interactions are defined as fol- 



lows: {eu,egg,eig) = (3,1,^0^^ 



1.73) kJ/mol, 



i<yii,crgg,aig) = (0.34,0.5,(CTn + crsg)/2 = 0.42) nm. 
The simulations are carried out in the NPT-ensemble 
(constant number of particles, pressure, and temper- 
ature). A Berendsen pressure scaling algorithm was 
applied, and the temperature was kept constant using 
a velocity rescaling thermostating procedure [21]. 
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Tabic 1: Simulation parameters of the different simulations. The Initial Conditions (IC) type refers to the 
configurations shown in Fig. 2. 



Exp. 


^[nm] 










IC type (Fig. 2) 


Stable? 


I 


15 


333 


12339 


2.70-10"^ 


3.0 


a 


yes 


II 


15 


342 


12330 


2.77-10-2 


3.1 


a 


yes 


III 


15 


432 


12240 


3.53-10-2 


3.9 


a 


yes 


IV 


15 


531 


12141 


4.44-10-2 


4.8 


a 


yes 


V 


30 


333 


52489 


0.63-10-2 


1.5 


a 


no 


VI 


30 


342 


52480 


0.65-10-2 


1.5 


a 


no 


VII 


30 


432 


52390 


0.82-10-2 


1.9 


a 


yes 


VIII 


30 


832 


51990 


1.60-10-2 


3.6 


a 


yes 


I-b 


15 


332 


12340 


2.69-10-2 


3.0 


b 


yes 


Ill-b 


15 


436 


12236 


3.56-10-2 


3.9 


b 


yes 


V-b 


30 


333 


52489 


0.63-10-2 


1.5 


b 


no nucleation 


Vll-b 


30 


432 


52390 


0.82-10-2 


1.9 


b 


no nucleation 


Vlll-b 


30 


837 


51985 


1.61-10-2 


3.7 


b 


no nucleation 




(Ng-333) gas atoms 
. (dissolved) ^_ 



Figure 2: The two types of initial conditions for the 
simulations, a) A preformed bubble containing 333 
gas-atoms surrounded by liquid. If there are more 
gas atoms in the system [Ng > 333) they are uni- 
formly dissolved throughout the liquid, b) All {Ng) 
gas atoms are uniformly dissolved throughout the liq- 
uid, so there is no prc-formed bubble. If the concen- 
tration of gas is high enough, homogeneous nucle- 
ation will occur forming a nanobubble. 



ration consists of a preformed bubble at a predefined 
radius i?o containing gas (333 atoms) and vapour. 
Outside the bubble the simulation box is completely 
filled with liquid, and the remainder of the gas is uni- 
formly dissolved throughout the liquid. For the sec- 
ond configuration, the simulation box is completely 
filled with liquid with the gas uniformly dissolved in 
this liquid (so no pre-existing bubble). In this con- 
figuration, a nanobubble will occur if the concentra- 
tion of gas in the liquid is high enough such that 
the energy barrier for homogeneous nucleation can be 
overcome. Since the pressure is maintained constant 
throughout the simulation, the box-size is allowed to 
vary to accommodate this. In practice, we find that 
the box dimensions never vary more than 10% from 
their initial values. The initial velocities for all atoms 
are sampled from a Maxwell-Boltzmann distribution 
at 300K. 



2.2 Initial conditions 



Results from the MD simula- 
tions 



We use two different atom start-position configura- 
tions, which are shown in Fig. 2. The first configu- 



A total of 8 different bubbles have been simulated (see 
Table 1) which started with a pre-existing bubble in 
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Figure 3: Snapshots from a selection of simulations. Note that the bubbles do not remain centered in the 
simulation box due to Brownian motion. This is not a problem since periodic boundary conditions are 
imposed, such that the bubble moves back into the simulation domain at the opposite site from which it 
leaves. Note also how the box size adapts to keep the pressure in the system constant. Left: Two simulations 
(III and VI) with initial conditions type a (Fig. 2). For simulation III, the bubble grows towards a stable 
radius. For simulation VI (larger £), the bubble completely dissolves within ~ 70 ns. Right: Bubbles 
simulated using initial conditions type b (Fig. 2). For simulation Ill-b, homogeneous nucleation occurs and 
the bubble grows towards the same equilibrium radius as the bubble in simulation III. For simulation Vlll-b, 
no nucleation occurs and the gas remains homogeneously dissolved throughout the liquid. 



the initial conditions (Fig. 2a). Of those configura- 
tions 5 additional simulations were performed using 
the initial conditions without a pre-existing bubble 
(Fig. 2b), to see whether the initial conditions affect 
the final result. The boundary of the bubble is de- 
fined at p* =0.5, where: 



Here, pv is the bulk number density of the gas/vapour 
phase inside the bubble and pi the number density in 
the bulk liquid. This boundary is then fitted with a 
circle giving R{t). Some snapshots of a selection of 
simulations are shown in Fig. 3. R(t) against time 
is plotted in Fig. 4 where we see that some bub- 
bles are stable, while others are not. As one would 
intuitively expect, the bubbles that are closest to- 
gether {£ = 15nm, configurations TIV) are stable, 
whereas some bubbles that are spaced further apart 
(£ = 30nm) are not. The stable bubbles benefit from 



their nearest neighbours, as they 'shield' the diffu- 
sive outflux that would normally lead to dissolution 
within microseconds. To confirm that these bubbles 
are truly stable, we extended one simulation (config- 
uration III) until t = 0.8 ^s, where we found that 
after i = 3 ns the radius R remained perfectly con- 
stant (see inset Fig. 4). 

The bubbles that are spaced farther from each 
other (£ = 30nm, configurations V-VIII) arc not 
always stable. The configurations with the least 
amount of gas (configurations V and VI) dissolve 
within 70 ns, whereas configurations VII and VIII 
- which contain more gas - are stable. 

What happens when the initial conditions are 
changed from a pre-existing bubble (Fig. 2a) to uni- 
formly dissolved gas in bulk liquid (Fig. 2b)? If a bub- 
ble forms, there is no reason why it shouldn't grow to 
the same equilibrium size as the corresponding bub- 
bles with different initial conditions. In fig. 5 we show 
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Figure 4: Results from the simulations with initial condition type a (Fig. 2). All bubbles in the small box 
{£ = 15nm) are stable (I- IV), whereas the simulations in the large box {£ = 30nm) and with little gas 
initially dissolved in the liquid (V,VI) are unstable. When a sufficient amount of gas is initially dissolved in 
the liquid (VII, VIII) the bubbles are stable even in the large box. The inset shows the extended simulation 
of configuration III, to verify that the bubble is indeed stable at long timescales. 



the results of some simulations with the alternative 
initial conditions (I-b, Ill-b, V-b, Vll-b, and Vlll-b) 
compared to the data of the similar bubbles with the 
original initial conditions (I, III, V, VII, and VIII). 
We see that when nanobubbles form, they indeed 
grow towards the same equilibrium size which they 
also achieve when starting with a finite size bubble. 
The reason that for these initial conditions in some 
cases no bubbles form is that the gas concentration 
(j) in the liquid is not large enough to overcome the 
nucleation barrier. 



4 Continuum description 

In this section, wc will use continuum fiuid mechanics 
to explain and predict the behavior of nanobubble 
clusters. In particular, we will address their stability 
and calculate their equilibrium size. We will however 
first address the subject of homogeneous nucleation. 



which is relevant for the simulations where there was 
no pre-existing bubble in the initial conditions. 

4.1 Nucleation theory 

In the case where there is no pre-existing bubble in 
the initial conditions (Fig. 2b), we are dealing with 
homogeneous nucleation, since there are no seeds 
(such as contamination) available to start heteroge- 
neous nucleation [22]. In the case of homogeneous 
nucleation, the change in free energy of the system 
when a bubble of radius R forms is: 

AG = ^nR^Gy + 4:TtR^j . (3) 

Here, AG is the energy gain or loss for the system 
to form a bubble of radius R. G„ is the (volumet- 
ric) energy associated with a unit volume of gas and 
is a negative number, hence it promotes nucleation. 
The liquid-vapour surface tension 7 is always posi- 
tive and therefore acts agains nucleation. For small 
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Figure 5: Bubble radius evolution as a function of 
time comparing similar systems with different initial 
conditions. When bubbles nucleate, they grow to- 
wards the same equilibrium size independent of ini- 
tial conditions. Although a stable bubble exists for 
configurations VII and VIII, the bubbles in configu- 
rations Vll-b and Vlll-b do not reach this state as 
the bubbles do not nucleate. Naturally, nucleation 
also does not occur for the configuration where no 
stable bubble can exist (V, V-b). 



4.2 Growth dynamics 
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Figure 6: Bubble radius function of time dur- 
ing growth. The data shown represents bubble VIII, 
as it grows towards its equilibrium size and has been 
shifted such that R{t = 0) = 0. At intermediate 
times an exponent of 0.56 is observed, consistent with 
diffusive bubble growth. In the final stage the avail- 
ability of gas is too low to sustain the growth rate, 
and the bubble settles at its equilibrium radius. 



R, the surface energy term usually wins, but when 
the bubble reaches a critical radius R* = —2"f/Gy 
then dG / AR\fi^ii* < 0, meaning that the bubble will 
grow. In our case, we start out with a bubble of zero 
radius, but due to thermal fluctuations small bub- 
bles appear randomly throughout the system. By 
increasing the gas concentration in the liquid, the 
magnitude of Gy increases as it becomes more and 
more favourable for the system to have gas atoms in 
the gas phase than in the dissolved state. This de- 
creases the value of R* , until it is small enough that 
the spontaneously forming tiny bubbles are already 
large enough {R > R*) to overcome the surface en- 
ergy penalty and grow. The growth stops when equi- 
librium between the gas and liquid phase is achieved 
[Eq. (5)]. 



The stage between nucleation of a gas bubble and it 
reaching its final size is governed by diffusive bubble 
growth. The relation between R and t for diffusive 
bubble growth is well known, namely a square root 
power law [23, 24], 

Rit)^iDt)K (4) 

where D is the diffusion constant of the gas in the 
liquid. In Fig. 6 we show the growth of bubble VIII 
on a log-log scale. The curve is shifted in time such 
that R{t = 0) = 0. A power law behaviour with 
exponent 0.56 is observed. This is very close to the 
expected exponent of i, suggesting that the growth 
of the gas bubbles in the simulation is indeed limited 
by diffusion. The power law of Eq. (4) is derived as- 
suming that the magnitude of the Laplace pressure is 
small compared to atmospheric pressure. This con- 
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dition is not fulfilled in particular for the tiny bub- 
bles in the beginning, which presumably accounts for 
the observed deviation from Eq. (4). Also, the final 
stage of bubble growth (t > 10 ns) exhibits a different 
behaviour with time, which is caused by the closed 
nature of the simulation system and thus the limited 
amount of gas available. In this stage the available 
gas in the system is depleted, and the bubble assumes 
its (final) equilibrium radius. 

4.3 Equilibrium radius 

Finally, the bubble reaches its equilibrium radius. 
When a bubble with radius R exists inside an infinite 
body of liquid, the concentration of dissolved gas just 
outside the bubble (r = R) is given by Henry's law: 



P - h 

-f — f^H 



(5) 



Here, p is the partial pressure of a specific gas in the 
gas phase, cf) the gas concentration of that specific gas 
inside the liquid, and kjj is Henry's constant. In this 
work we will use the (dimensionless) mole fraction of 
the gas and liquid as the concentration (p. The total 
pressure in a (2D) bubble is given by: 



Pb 



P' + R 



(6) 



where pi is the pressure in the surrounding liquid, 
which is usually negligible in the case of nanobubbles. 
From this relation, we can see that surface tension is 
a strong driving force for dissolution, especially as 
the bubble radius R becomes small. 

kn is Henry's constant and depends on the type 
of liquid and the type of gas involved. kn is 
temperature-dependent, but in this work we will con- 
sider a system with a fixed temperature. In sepa- 
rate measurements (in a system consisting of a liquid 
phase in equilibrium with a gas phase at p = 1 atm) 
we found that kn ^ 10^ Pa for the Lennard- Jones 
fluids considered in this work. 

A gas concentration gradient induces a diffusive 
mass flux J according to Fick's law: 



In the case of a cluster of bubbles, there is a limited 
amount of liquid present between neighbouring bub- 
bles. Since the amount is limited (and small if the 
bubbles arc sufficiently close to each other) the con- 
centration is noticeably affected by the gas flow out of 
the bubbles. Eventually, as more and more gas enters 
this space, the gas concentration at the mid-point be- 
tween these bubbles will reach the gas concentration 
at R as prescribed by Henry's law [Eq. (5)]. When 
this has happened, there is no concentration gradient 
any more (the concentration is equal everywhere) and 
no diffusive gas flow exists: the bubbles are stable. 

Using Henry's law we can predict the equilibrium 
size of these nanobubbles. For this, we have to ex- 
press the concentration of the gas in the liquid as a 
function of the bubble radius R. First, we know that 
the amount of gas in the liquid is simply given by: 



n'{R)^Ng-n'iR) , 



(8) 



where is the number of gas molecules dissolved in 
the liquid, Ng the total number of gas atoms in the 
system, and the amount of gas molecules inside 
the bubble. The amount of gas inside the bubble can 
be related to the size of the bubble, using the Laplace 
pressure and the ideal gas law: 



b , 27r7M 



(9) 



Here, k^ is Boltzmann's constant, T the temperature 
of the system. To obtain this relation, wc assumed 
that the Laplace pressure difference is much greater 
than the ambient pressure, which is a valid assump- 
tion for nanobubbles exposed to atmospheric ambient 
pressure. Henry's law [Eq. (5)] dictates the equilib- 
rium condition, so by combining Eqs. (8) and (9) we 
can now solve for the equilibrium bubble radius R^q. 



4:7r^d 



1± Wl-4 



kBTknN^ 



(10) 



(7) 



Since the equation is quadratic we get two solutions, 
of which we cannot say a priori which one is valid, 
because both arc positive and finite. 
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Figure 7: Critical number of gas atoms N* against 
Ni. The critical number of gas atoms is the amount 
of gas atoms required to prevent full dissolution of 
a nanobubble. As the system is larger (hence larger 
separation distance between nanobubbles) more gas 
is required to sustain a stable nanobubble. 



4.4 Comparison with simulations 

We can now compare the simulation results to the 
model [Eq. (10)]. 

From Eq. (10) it is clear that for a physical solution 
to exist the number of gas atoms Ng must be larger 
than: 



ksTk 



H 



(11) 



Hence, N* is a critical amount of gas atoms: When 
fewer gas atoms are present in the system the bubble 
will completely dissolve. 

In Fig. 7 this critical gas content as a function of 
Ni , which is proportional to the system size is shown. 
In the figure the two system sizes considered in this 
study are annotated. From the figure it is clearly that 
for larger systems (hence, larger nanobubble sepera- 
tions and higher Ni) the required amount of gas in- 
creases. Intuitively, one would expect that the ratio 
N*/Ni is constant, but from Eq. (11) and Fig. 7 we 
can see that this is not the case. Instead, the mini- 



mum number of gas molecules scales with the square 
root of the amount of liquid present, meaning that it 
scales linearly with the interbubble distance £. 

The predicted equilibrium radius [Eq. (10)] as a 
function of gas content for the small systems {£ = 
15nm) and the large systems {i ~ 30 nm) are shown 
in Fig. 8(a). The measured equilibrium radii are 
also depicted, and show excellent agreement with the 
model. For configurations V and VI (large systems), 
the model correctly predicts that there is no stable 
bubble size, hence they fully dissolve. For configu- 
rations Vll-b and Vlll-b the model predicts that an 
equilibrium radius does exist, however in these cases 
there are no bubbles since the nucleation barrier is 
too high, which is not accounted for in the model. 
By nondimensionalizing the equilibrium radius R^q 
with £, and rescaling Ng with N* we can collapse the 
data on a single curve, independent of system size, 
see Fig. 8(b). 

For Ng :$> N*, Eq. (10) can be expanded as: 




(12) 



These asymptotic solutions are plotted in Fig. 8b as 
the black dashed lines. The first solution of Eq. (12) 
represents the limit where all gas atoms are contained 
within the bubble. This can also be seen by compar- 
ing this result with Eq. (9). It is apparent from Fig. 8 
that almost all bubbles are very close to this limit i.e. 
there is virtually no gas dissolved inside the liquids 
in the simulations. 



5 Nanobubble cluster stability 

Finally, we will discuss the stability of the entire clus- 
ter of nanobubbles. Although it is now clear that 
bubbles that are surrounded by mirror images can 
indeed be in equilibrium, this does not mean that the 
entire cluster is stable. There are two obvious threats 
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Figure 8: (a) Equilibrium radii Rf,q as measured in the simulations (symbols) and predicted by the model 
(lines) as a function of Ng. The blue data correspond to the small system {£ = 15nm), whereas the red 
data correspond to the large system {£ = 30nm). The model nicely predicts the equilibrium size in both 
systems, but of course does not account for the configurations where no bubble nucleation occurs (V-b, VII- 
b, and Vlll-b). To obtain the theoretical curves, a value of kn = 1-5 • 10^ Pa was used as a fit parameter. 
Independent measurements have shown kn ~ 10^ Pa. (6) Collapse of the nondimensionalized equilibrium 
radii. All data points nicely follow the master curve. The simulations where no nucleation occurs arc not 
included in this plot, to emphasize the good prediction of both i?eg and TV* by the model. The dashed lines 
indicate the expansion when iV^ S> N*, Eq. (12). 



to this stability: (i) the Brownian motion of bubbles 
can cause neighboring bubbles to collide, leading to 
coalescence and coarsening of bubbles within the clus- 
ter and {ii) when neighboring bubbles are not exactly 
equally sized, smaller bubbles would drain into larger 
bubbles via diffusion (similar to Ostwald ripening). 

First, we address the possibility of bubbles collid- 
ing. Since the distance between neighbouring bub- 
bles is at best ten times as large as the bubble ra- 
dius. Brownian motion of nanobubbles can indeed 
lead to collisions (see Fig. 3, where the nanobubbles 
are shown to move around). An obvious way to pre- 
vent this from happening is to make sure that the 
bubbles repel each other. This repulsion could be 
electrostatic (e.g. by using an ionic surfactant, or by 
the intrinsic negative charge of air bubbles in water 
[25-27]), which is sufficiently long-ranged. Adding 



salt to the solution, which leads to screening of the 
electrostatic fields around the nanobubbles, would 
then reduce the stability of bulk nanobubble clusters 
which has indeed been observed in experiments [3]. 



The second issue (Ostwald ripening) is harder to 
prevent, and we expect that in time coarsening would 
indeed occur. Of course, the larger bubbles that are 
formed in this way are still able to provide some 
shielding, but the polydispersity of bubble sizes will 
make a theoretical analysis difficult. Therefore, at 
this point, we do not have an explanation for the 
measured stability by Ohgaki et at [2] where Ost- 
wald ripening apparently has been suppressed. 
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6 Conclusions 

In conclusion, we have shown that bulk nanobubble 
clusters can indeed be stable under specific condi- 
tions. First, the individual bubbles are surrounded 
by similar nanobubbles. The distance between these 
bubbles must be small enough, such that the bub- 
bles can succeed in saturating the liquid between the 
bubbles with gas, before the bubbles are completely 
drained. When a cluster of stable bubbles exist, and 
the system is closed (so no gas can escape) the bub- 
bles will in principle live forever, as long as bubbles 
cannot merge and Ostwald ripening is somehow pre- 
vented. The merging of bubbles can be prevented by 
ionic surfactants. Ostwald ripening is harder to pre- 
vent in theory, but experimental results showing that 
nanobubble cluster can indeed be stable for longer 
times indicate that there exists a mechanism that can 
prevent Ostwald ripening from occuring [2-4] . 
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